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1 Abstract 

Particle shape and aggregation have a strong influence on the spectral profiles of infrared phonon 
bands of solid dust grains. In this paper, we use a discrete dipole approximation (DDSCAT; 

a cluster-of-spheres code following the Gerardy-Ausloos approach (MQAGGR; [2]) and a 
T-matrix method (SCSMTM; for calculating IR extinction spectra of aggregates of spherical 
silicon carbide (SiC) particles. We compare the results obtained with the three different methods 
and discuss differences in the band profiles. 



2 Introduction 

Grain growth by aggregation is an important process in dense cosmic environments as well as in 
the earth's atmosphere. Besides influencing dynamic properties it also changes the absorption 
and scattering properties of the solid dust particles for electromagnetic radiation (e.g. [1]). This 
is especially true in spectral regions where resonant absorption occurs, such as the phonon bands 
in the infrared. It is quite well known that shape and aggregation effects actually determine 
the band profiles of such absorption and emission bands, which hinders e.g. the identification of 
particulate materials by their IR bands, but detailed investigations especially of the infiuence of 
grain aggregation are still lacking. We plan to set up a spectroscopic experiment (see Tamanai 
et al., this volume) for measuring aggregation effects on IR extinction by dust particles dis- 
persed in air and, simultaneously, have started to use light scattering theory in order to predict 
numerically the band profiles for different aggregation states. 




Figure 1: The clusters frac7, lin9 and sc8. 
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Figure 2: The refractive index for /5-SiC with the damping constant 7 = 2 cm~ . The damping 
constant 7 (inverse proportional to the phonon hfe time) is an "ad hoc" introduced parameter, 
which in a perfect crystal reflects the anharmonicity of the potential curve. High-quality SiC is 
characterized by damping constants of 1-3 cm~^. 

3 Structure of the clusters 

We consider three-dimensional clusters of identical touching spherical particles arranged in 
three different geometries: fractal, cubic, and linear. For a high precision of the calculations, 
we restrict the number of particles per cluster to less than 10. Therefore, we have selected 
only three geometries, namely a "snowflake I'st-order prefractal" cluster (fractal dimension 
D = In 7/ In 3 = 1.77, where one sphere is surrounded by six others along the positive and 
negative cartesian axes, a cluster of eight spheres arranged as a cube and a linear chain of nine 
spheres. All of the clusters (Fig. 1) consist of spheres with radii R = 10 nm and are embedded 
in vacuum (or air). 

For the optical constants of the particle we have chosen the data of /3-SiC in the wavelength 
range 10-13^m, calculated from a Lorentzian oscillator- type dielectric function describing the 
phonon resonance in this wavelength range (see |S]). On the one hand, this phonon resonance is 
of practical importance since it is observed as an emission band from dust particles in carbon star 
envelopes. On the other hand, depending on the resonance damping parameter, it represents 
a model material of a very high complex refractive index with |m| > 10 and sharp surface 
resonances in the wavelength range between the LO and TO frequencies. 

4 Computational approaches 
4.1 The DDA method 

The discrete dipole approximation (DDA) method is one of several discretisation methods (e.g. 
[Zj, |H]) for solving scattering problems in the presence of a target with arbitrary geometry. In 
this work we use the DDSCAT code version 6.1 _1,, which is very popular among astrophysicists. 
In DDSCAT the considered grain/cluster is replaced by a cubic array of point dipoles of certain 
polarizabilities [S]. The cubic array has numerical advantages because the conjugate gradient 
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Figure 3: Extinction profiles in the SiC phonon band for optical constants with 7 = 2 cm^ and 
two different cluster geometries, the linear chain of 9 touching spheres and the cubic cluster 
of 8 touching spheres. Dotted curves show the results obtained with DDSCAT, dashed curves 
with SCSMTM, and sohd curves with MQAGGR. 

method can be efficiently applied to solve the matrix equation describing the dipole interactions. 
By specifying an appropriate grid resolution, calculations of the scattering and absorption of 
light by inhomogeneous media such as particle aggregates can be carried out to in principle 
whatever accuracy is required. For the sc8 cluster we used a grid of 36 x 36 x 36 dipole, which 
provides 23 752 dipoles in the cluster, while for the lin9 cluster a grid of 64 x 64 x 64 dipole 
was used providing 1676 dipoles in the cluster. The program performs orientational averaging 
of the clusters. 

4.2 The clusters-of-spheres method 

The clusters-of-spheres calculations have been performed using (1) the program developed by M. 
Quinten (MQAGGR, commercially available) based on the theoretical approach by ^H] and (2) 
using the T-matrix code by D.W. Mackowski (SCSMTM) calculating the random-orientation 
scattering matrix for an ensemble of spheres |3j. Both programs aim at solving the scattering 
problem in an exact way by treating the superposition of incident and all scattered fields, 
developed into a series of vector spherical harmonics. Available computer power, however, forces 
to truncate the series at a certain maximum multipole order npolmax, which in both programs 
can be specified explicitly. Furthermore, both programs perform an orientational average of 
the cluster, the resolution of which was set to 15 degrees in MQAGGR and 10 degrees in 
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SCSMTM for theta (the scattering angle). The variation in the azimuthal angle is not specified 
in SCSMTM. In MQAGGR it is varied betwen and 360 degrees, again with a resolution of 15 
degrees. 

5 Results and discussion 

The extinction profiles shown in Fig. 3 display the results obtained for two cluster geometries 
and three different codes. The profiles obtained with the two clusters-of-spheres codes at the 
same multipolar order (13 for the linear chain and 15 for the cubic cluster) give quasi-identical 
results in the case of the linear chain but differences for the cubic cluster. SCSMTM gave 
less resonances than MQAGGR, which is true for all lower maximum multipolar orders as 
well. The reason for that could possibly be different ways of orientational averaging. This is not 
completely understood yet. Generally speaking, both codes didn't converge up to the multipolar 
orders tried (up to 2h of CPU time per wavelength point). 

DDSCAT at the maximum resolution used (1676 dipoles in the linear cluster, 23752 dipoles 
in the cubic one) gave much smoother profiles, i.e. with less distinct resonances. DDSCAT 
possesses limitations for large refractive indices (|llj). On the other hand, the resonances pro- 
duced by the clusters-of-spheres codes depend so much on the maximum multipolar order that 
it seems reasonable to assume that the (unreached) converged spectrum would show a more 
smooth profile as well. Results with lower dipole resolution tend to show sharp resonances as 
well. 
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